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Abstract 

We introduce a model of binary data corruption induced by a Brownian agent 
(active random walker) on a ci-dimensional lattice. A continuum formulation 
allows the exact calculation of several quantities related to the density of cor- 
rupted bits p ; for example the mean of p, and the density-density correlation 
function. Excellent agreement is found with the results from numerical sim- 
ulations. We also calculate the probability distribution of p in d = 1, which 
is found to be log-normal, indicating that the system is governed by extreme 
fluctuations. 

PACS numbers: 05.40.+j, 66.30. Jt, 82.30.Vy 
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I. INTRODUCTION 



Brownian motion is one of the fundamental processes in Nature. Originally observed in 
the irregular motion of pollen grains by the botanist Brown Q, and cast into the language 
of the diffusion equation by Einstein [fj], it has now been applied, in the mathematical 
framework of random walks 0, to an enormous variety of processes in the physical sciences 
and beyond. A very rich field of research has been built up around the behavior of a random 
walk coupled to a disordered environment , a good example being the anomalous diffusion 
of electrons in a disordered medium |]5| . Also, one can consider a random walker being the 
disordering agent in its environment. Applications of the latter include the tagged diffusion 
of atoms in a crystal or magnetic disordering due to a wandering vacancy §|. 

In this paper we will introduce a particularly simple example of an active random walker 
(or Brownian agent) disordering its environment. Although the model is interesting in its 
own right, we believe it will have a useful application to the study of data corruption in 
ultra-small storage devices. Before pursuing this connection, we shall briefly describe the 
model (which will be more carefully defined in the next section). The two main features 
of the model are first, the Brownian agent (BA) performs a pure random walk - it is not 
affected by the environment in any way. Second, the environment is bistable. That is to 
say, it is composed of elements which may only exist in one of two possible states (see Ref. 
ToU for a loosely related random walk process). Thus we can consider the environment to 



be composed of binary data (our favored realization), magnetic spins, chemical species A 
and B, and so on. As the BA wanders through the environment it has a certain probability 
to switch the value of an element in its immediate vicinity. Thus, if we start with a system 
in which all elements exist in the same state ('up' say), and introduce the BA at the origin, 
then after some time, there will be a region around the origin in which the elements will 
be found in a mixture of 'up' and 'down' states. Naturally, the linear size of the region 
will grow on average as y/i. The more subtle question concerns the degree of disordering 
which exists for elements within this region, and also their spatial correlations.. As we 
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shall see, the statistics of the disordered elements are very rich. This is most convincingly 
demonstrated by the dominance of extreme fluctuations; for instance, the distribution of 
disordered elements is log-normal. Thus, typical and average events are quite distinct, and 
become ever more so as time proceeds. 

Before giving an outline of the paper we shall say a few words about the potential rel- 
evance of this system to data corruption. With the advent of semiconductor memories 
(for dynamic random access memories (DRAM's) and various types of read only memories 
(ROM's)), there has been a tremendous drive within the semiconductor industry to produce 



ever-smaller memory devices |IT|,|12|. There are many properties (e.g. stability, power con- 
sumption, volatility, and cost) which must be balanced in the design of such devices. These 
factors determine the type of material used, and the geometry, dimension and architecture 
of the device. (For instance, three dimensional arrays have a very efficient address structure, 
and are stable against interference from bombarding alpha particles, but are very expensive 



to produce ||11|| .) One of the main issues is the stability, or reliability, of the device. In 
semiconductor memories, there are many physical effects which can create hard errors (de- 
struction or corruption of the device itself) or soft errors (corruption of the data stored in the 
device). In the latter category, the most common problems originate from electron clouds 
caused by alpha particles, but soft errors may also arise from electromigration and charge 



diffusion [R|. The key point is that different corruption mechanisms operate on different 



time scales (leading to the famous bath-tub curve of device reliability |T2[)- It is therefore 
important to know on what time-scales one should expect significant corruption from a given 
process. The model we propose here (namely data corruption via a BA) is probably not 
relevant for today's semiconductor devices, since there are so many 'mesoscopic' processes 
occurring on the level of a flip-flop, that subtle correlations due to a BA will be washed out. 
However, we can look ahead to the new generation of (quantum) storage devices, in which a 
single electron (controlled in a gate via coulomb blockade) can store one bit of data. In this 
case, a microscopic BA may indeed play an important role in data corruption, and it will be 
necessary to understand its time-scales and efficacy of operation, so that we can minimize 



its influence. This paper constitutes a first step towards gaining such an understanding. 

The outline of this paper is as follows. In the next section we shall carefully define the 
model (using discrete space and time) within the master equation formulation of stochastic 
processes. We shall derive some general statistical properties of the process, but we shall 
not enter into any explicit calculations. This is deferred to sections III-V in which we 
introduce a very simple continuum theory for the process, which is motivated by viewing the 
process as a stochastic cellular automaton. In section III we derive this continuum theory 
and, using the complementary descriptions of quantum mechanics (i.e. the Schrodinger 
equation and the Feynman path integral), we shall demonstrate its equivalence to the master 
equation formulation of section II. We then examine the case of one spatial dimension 
in section IV. The model is tractable, and a great deal of information may be derived 
concerning the mean density of disordered elements, their spatial correlations, and finally 
their entire distribution function. In section V, we briefly study higher dimensions, and 
derive some general statistical properties for the process, for an arbitrary spatial coupling 
between the BA and its environment. We also derive an expression for the mean density 
of disordered elements in two dimensions. In section VI we present results from extensive 
numerical simulations of the discrete process. In all cases, we find good agreement between 
the simulation results and the predictions of the continuum theory. We end the paper with 
a summary of the work, and some ideas for future study. 

II. DISCRETE FORMULATION OF THE MODEL 

We consider binary data bits on a <i-dimensional hypercubic lattice. For convenience we 
shall represent each bit by an Ising spin a r , where the index r represents a discrete lattice 
vector. The spin takes the value +1(— 1) for a data bit which is uncorrupted (corrupted). 
Thus, the initial configuration is a lattice of spins, all of which take the value +1. [We 
prefer to describe the system almost exclusively in terms of the spin variables. Thus we 
shall use phrases such as 'magnetization density', or 'global magnetization'. The translation 
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of these quantities to the corresponding properties for corrupted data bits is immediate as 
one only need replace a r by (1 — 2n r ), where n r denotes the presence (with value unity) or 
absence (with value zero) of a corrupted bit. Similarly, we shall often refer to the average 
magnetization density m , which is related to the average density of corrupted bits p via 
m = 1 — 2p .] We denote the position of the BA by the lattice vector R. Each time step, 
the BA has a probability p to make a jump to one of its (2c?) nearest neighbors. For the 
sake of generality, we will not insist that the BA always flips a spin (i.e. changes a data 
bit) as it moves. Thus, on a given jump, we allow the BA to flip the spin at the site it 
is leaving, with a probability q. In this section we shall describe the process via a master 
equation |T|| . Namely, we shall define the dynamics through the evolution of the distribution 
P(R, {er r }, t), which is the probability that at time t the BA is at position R, and the spins 
have configuration -R}. Given the above rules, the master equation takes the form 

P(R, R}, t + 5t) = (l- p)P(R, R}, t) + E P(R + 1, R}, *) 

+ §j£P(R + l,---,-<TR + i,---,t) , (1) 

where {1} represent the 2c? orthogonal lattice vectors (which have magnitude /). 

The natural quantities to extract from the distribution P are conditional averages. The 
simplest is the mean value of the spin at position ri at time t, given the BA is at position 
R. This is defined via 

0(ri, t I R) = Tr CT a ri P(R, R}, t) . (2) 

Higher order conditional averages may be defined accordingly. Performing the spin trace 
over the master equation with a weight of a r yields 

Q(r,t + 5t | R) - 0(r, t | R) = Y d E[ Q ( r ' * I R + 1) - 6(r, t | R)] 

-f e(r,t|r)]p riR+1 . (3) 

The above equation has a physically appealing form. The rate of change of G has two 
contributions. The first is lattice diffusion, as given by the first sum on the r.h.s. The 
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second contribution vanishes unless the spin in question is in the immediate vicinity of the 
BA, in which case it acts as a sink. 

At this point in the discussion it is worthwhile to consider the continuum limit. Namely 
we take the time scale 5t and the lattice scale I to zero, and define a diffusion constant 
D = l 2 p/5t. We also introduce a coupling A oc ql/5t. Then replacing the Kronecker 5- 
function in Eq.(H) by a Dirac 5- function we find 

<9 t 6(r, t | R) = DV 2 R <d(r, t | R) - A 6(r, t | r)5 d (r - R) . (4) 

It is important to note that this continuum equation is not strictly derived from Eq.(||), as 
we have not proved that the continuum limit exists. In fact, we shall find that for d > 2, 
the lattice scale / is crucial, and consequently we must soften the Dirac ^-function to a 
function A;(r) which is sharply peaked (over a linear scale I) with unit integral over lZ d . One 
actually expects this to be the Eq.(|j) is the imaginary-time Schrodinger equation 

for a particle under the influence of a repulsive 5-function potential. [Note the independent 
spatial variable in this quantum system is R, with the variable r simply labelling the position 
of the potential.] It is well known [JTSJ that the repulsive ^-function potential is 'invisible' to 
the particle for d > 2, and one usually cures this by smearing the potential just as described 
above. The quantum mechanics analogy will prove useful in the next section when we 
construct an alternative continuum model. 

Before leaving this section we shall indicate the derivation of a non-trivial statistical 
relation hidden inside Eq.([3]). First, we must define the initial condition. As mentioned 
before, given that up-spins denote uncorrupted data, the initial value of each spin is +1. 
There is a slight subtlety of definition regarding the value of the spin at the site where the 
BA is initially planted. (This position shall be taken to be the origin, without any loss of 
generality.) We shall take this spin to be initially —1 so that immediately after the BA has 
moved away the spin at the origin has value +1. Thus we have 

P(R,K},0) = Vo^ o ,-iII^,i , (5) 

r^O 
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and consequently, 

e(r,0 I R) = 5 R ,o(l " 25 r ,o) • (6) 

We may obtain the following two averages from 0(r, t | R). The first is the average 
value of the spin at the origin. This is simply given by m (i) = J2n 0(0, t | R). The second 
average is the quantity Xr ©(R> t | R), which corresponds to averaging the value of the spin 
at the site where the BA happens to be at time t. One can prove that 

£e(o,t | R) = £e(R,t | R) , (7) 

R R 

for all t. We arrive at the above result by essentially solving the partial difference equation 
(0) using discrete Fourier and Laplace transforms. The details can be found in Appendix A. 

This result is useful for proving a more physically relevant relation. Let us denote the 
average global magnetization by 

M W = EE[0MIR)-0MIR)]> (8) 

R r 

where we have defined it relative to the (infinite) initial magnetization. This quantity 
essentially measures the average of the total number of corrupted bits (up to a factor of 2). 
Summing Eq. (Q) over r and R gives 

M(t + 5t) - M(t) = 2pq 9(r, t | r) . (9) 

r 

Then using (0) we can rewrite the above relation in the form 

M(t + 5t)-M(t) = 2pqJ2®(0,t | R) = 2pq m (t) . (10) 

R 

In other words, the rate of change of the mean global magnetization is proportional to the 
mean magnetization density at the origin. This is a non-trivial relation between a global 
and a local quantity. 

In principle, one can obtain exact results for many interesting quantities (like the mean 
magnetization density, or correlation functions) by directly solving for the conditional aver- 
ages, as illustrated in Appendix A. However, we prefer to obtain results from a continuum 
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theory; partly because the calculations are a little easier, but more importantly because we 
can access more sophisticated properties of the system, such as the probability distribution 
of the coarse-grained magnetization density. 

III. CONTINUUM THEORY 

In this section we shall motivate a particularly simple continuum description of the data 
corruption process, and show its equivalence to the discrete theory of the previous section. 

There is an alternative method of characterizing the evolution of the system, other than 
using the evolution of the probability distribution P(R, {a r },t) via the master equation. 
This method consists of writing the local rules for the process in the spirit of a stochastic 
cellular automaton (SCA) [[Uj. Let us focus on the case that at each time step the BA 
makes a random jump to one of its nearest neighbors, and that the spin at the site which 
it leaves behind, definitely flips. This corresponds to setting p — q — 1. The local rules for 
such a process are easily written down. Let us denote the time-dependent position of the 
BA by R(t), a randomly chosen unit lattice vector by and the time-dependent value of 
the spin at site r by <J r (t). Then we have 

R(t + 5t) = R(t) + l(t) (11) 
a r (t + 5t) =a r (t)(l-25 rm) ) (12) 

We are interested in a continuum limit of these two rules. The first is nothing more 
than a random walk. We take the lattice vector R(£) to be a continuum vector (i.e. each 
of the d components is a real number), and we replace the random unit lattice vector l(t) 
by a continuum vector £(t), each component of which is a uncorrelated Gaussian random 
variable with zero mean (i.e. ^(t) is a white noise process). The correlator of £ is given by 

where here and henceforth, angled brackets indicate an average over the noise (or equivalently 
the paths of the BA). Then, on taking 8t — > 0, Eq.flTTl) assumes the form 
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£ = ««>, (14) 
which is the familiar equation for a continuum random walker where D' is the diffusion 



constant |Tj]. The second SCA rule is more complicated to generalize to the continuum. As 
a first step let us define a coarse-grained magnetization density <p(r,t) in the following way. 
We imagine defining a large region around the lattice site r and summing all the spins in 
that region. Their sum (suitably normalized) constitutes 0(r, t), with the label r denoting 
a point in the lZ d continuum. An entirely analogous procedure is used in motivating the 



Landau- Wilson free energy functional from the Ising model of ferromagnetism [I7|. The 



difficulty in our case is that we cannot derive a closed equation for from the discrete rule 
(|l2l). We therefore make the following approximation. Splitting the r.h.s. of (0) into two 
pieces, we see that the first may be taken over to the l.h.s. which may then be taken to be 
a time derivative in the limit of 5t — > . The second piece resembles a decay term centered 
at r = R. So, we postulate that the coarse-grained magnetization density satisfies 

dt<f>{T } t) = -A>(r, t)6 d (r - R(t)) , (15) 

where A' is a phenomenological parameter which describes how strongly the magnetization 
density is coupled to the BA. We stress that the field <p(r, t) is a function of the continuous 
space and time variables r and t, and a functional of the path R(£) of the BA. 

Now, the above heuristic derivation of the continuum theory was based on a SCA for the 
case in which the BA always moves (p = 1), and for which the spin located at the previous 
BA position is always flipped (q = 1). In general p and q are both less than unity. Intuitively 
we expect a very simple renormalization of our phenomenological parameters as p and q are 
changed. The diffusion constant D' should be proportional to p, and the strength of the 
spin-BA coupling A' should be proportional to both p and, more importantly, q. Thus we 
see a very close correspondence between D' and A' in the current continuum theory, and 
the parameters D and A which were introduced in the continuum limit of the discrete 
equation (|3|). In fact they are identical, as will emerge in the following discussion. 
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One of the positive features of the continuum theory as described by ([15]) is that one 
may immediately integrate the equation to find the magnetization density as an explicit 
functional of the path of the BA. As an initial condition we take 0(r, 0) = 1 for all r. The 
subtlety encountered in the discrete theory concerning the initial value of the spin at the 
origin disappears here since the coarse grained function is not sensitive to the value of one 
inverted spin. Straightforward integration of (|i~5|) yields 



It is important to note at this stage that the magnetization density is clearly positive 
for all r and t. Therefore within our continuum formulation, we have ignored paths of the 
BA which create large patches containing a majority of down spins {i.e. corrupted bits). 
Such patches will occur, but their frequency of occurrence is certainly very small since the 
system starts in a completely uncorrupted state. For instance, the probability for the BA 
to create a purely negative domain of iV spins is of the order e~ N . Therefore, so long as we 
coarse-grain the original spin model over a sufficiently large scale, we can be confident that 
the most important configurations have been retained in the continuum theory. Ultimately, 
one must justify such an approximation a posteriori by comparison with either results from 
the discrete theory, or from numerical simulations. As we shall see, both of these support 
the current continuum model and the approximations contained therein. 

We shall now connect the continuum theory as given by (|T6|), and the continuum limit 
(|j) of the discrete theory, as given by (|J). The mean magnetization density in the discrete 
theory is given by 



where in the continuum limit, we have replaced the sum over BA positions by an integral, 



alternative continuum theory, we can find the mean magnetization density by averaging the 
coarse-grained density 0(r, t) over all paths R(t). Each path is weighted by a Gaussian 
factor 




(16) 




R 



(17) 



and the field O satisfies the imaginary-time Schrodinger equation as given in (|]). In the 
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^{-^/tf { (0'}=^{-^/* (§)] , (18) 

where M is a normalization factor. Therefore we can write the average of as a functional 
integral 

m(v,t) = <0(r,t)> =UJ VR(t') exp j-^ j dt' (^j \ <j>(r,t) , 

= MjvR(t')e W ^-jdt' ^{^J ' + \'5 d (r-W)) J , 



R, 



A/" / d d i2/ / DR(t') exp J - j dt' 

I o 



dR 



2 



+ A'5 d (r - R(t')) 



2D' V dt' 

where we have used ([Tj^) in going from the first line to the second, and we have introduced 
the final position of the BA {i.e. R(t)) as a free integration variable Rj in rewriting the 
second line as the third. The reason for this cosmetic change is to make explicit the fact that 
(0(r, t)) can be expressed as a spatial integral over the final BA position, where the integrand 
is itself a path integral over BA trajectories. This path integral is nothing more than a re- 
expression of the solution of an imaginary time Schrodinger equation (using the well-known 
Feynman path integral formulation of quantum mechanics for a particle in a repulsive 
5-function potential. We can now see the connection: Eq. flT9| ) is an exact restatement of 
Eqs. (f§) and (0), with the identification D' = D and A' = A. [So, henceforth we shall 
drop the primes in the material parameters.] To summarize, by utilizing the complementary 
formulations of quantum mechanics via the Schrodinger equation and the Feynman path 
integral, we have shown that the continuum limit of the master equation is identical to the 
continuum theory constructed at the beginning of this section. 

Two final points are in order. First, as noted in the previous section, the ^-function 
potential must be replaced by a smeared function Aj(r) for d > 2; thus in our continuum 
theory encapsulated in Eq. (|T5D , we shall make a similar replacement when studying two or 
higher dimensions. Second, we have established a connection between the master equation 
and Eq.(P^[) only at the level of the first moment. It is straightforward to extend each formu- 
lation to higher order correlation functions, and indeed one finds an exact correspondence. 
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, (19) 



For instance, we can define the conditional spin-spin correlation function within the discrete 
theory 

e(ri, r 2 , t | R) = Tr CT a ri a r2 P(R, {a r }, t) . (20) 

Using the master equation one can show that in the continuum limit this function satisfies the 
Schrodinger equation for a particle under the influence of two repulsive ^-function potentials 
located at ri and r 2 . Similarly, we can construct the coarse-grained two-point correlation 
function from Eq.([L£j) by evaluating (0(ri, t)0(r 2 , £)). It is easy to see that this quantity is 
given by an integral over the analogous path-integral for two repulsive 5-function potentials. 

Having completed our formulation of a simple continuum theory, and shown its equiva- 
lence to the continuum limit of the master equation, we shall proceed to the next section in 
which we present a comprehensive solution of the model in one dimension. 

IV. RESULTS IN ONE DIMENSION 

In this section we restrict ourselves to one dimension. This does not necessarily mean 
a single chain of sites. Rather, we shall exclusively study the continuum theory of the last 
section, and in this case, for large enough times, d = 1 refers to any system which has an 
infinite longitudinal dimension, and finite transverse dimensions (for instance an infinitely 
long strip). This is the case, since as time proceeds, the correlation length will eventually 
become greater than the transverse size of the system, thereby only allowing the longitudinal 
fluctuations to continue growing, as is the case in a strictly one-dimensional system. 

The continuum model described in the previous section can be viewed as a 'non- 
conserved' cousin of the continuum theory of vacancy mediated diffusion (a process which 
in the spin language conserves magnetization) introduced recently ||. An exact analysis 
of the latter theory was possible using infinite order perturbation theory in the spin-BA 
coupling A. We shall use the same technique here, as it leads rather directly to a full so- 
lution. Alternatively, one may solve the Schrodinger equation for the conditional averages. 
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However, there are some important quantities (like the distribution of the magnetization 
density) which cannot be easily recovered from the latter approach. 

Our starting point is the integrated solution of the continuum formulation as given in 
Eq. (|i~6"|) . First, we shall derive an expression for the magnetization density m(x,t). Per- 
forming a direct average of Eq.flTBj) and expanding in powers of A (= A'), we have 



m(x,t) = (<f>(x,t)) = J2(- X ) H Xn(x,t) 
where Xo{ x i t) = 1, and for n > 



(21) 



n=0 



Xn(x, t) 



nl 



J dr 5(x - R(t)) 



(22) 



We refer the reader to Appendix B in which the above average is explicitly calculated. The 
result is 



C Tl >n— 1 

Xn(x, t) = J dn J dr 2 - ■ ■ J dr n g(0, n - r 2 ) • • ■ ^(0, r n _i - T n )g(x, r n ) 



(23) 



where g(x,t) = (2nDt)~ 1 ^ 2 exp(— x 2 /2Dt) is the probability density of the BA. 

The structure of Eq.(p3|) is that of an n-fold convolution, so we may utilize a Laplace 
transform to good effect. We have (for n > 0) 



Xn(x,s) = dte st Xn(x,t) = - g(0, s)" 1 g(x, s) . 



(24) 



Performing the sum over these functions as dictated by Eq . (pl|) we find 

\g(x, s) 



mix, s 



(25) 



1 + Xg(0, s) 

We note in passing that a similar result is easily derived for any d G (0,2). The case of 
d > 2 is more complicated as the function g(0,t) is no longer integrable. 

This expression for the Laplace transform of m is exact. This will prove to be important 
when we come to evaluate the distribution function of V(m). The inverse of the Laplace 
transform is given by 



m(x, t) = erf 



x 



{2Dty/ 2 



+ exp 



'Aid A 2 t\ 

1 erfc 

D 2D 



A 



\2D 



+ 



\x\ 



{2Dty/ 2 



(26) 
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where erf(z) and erfc(z) are error functions |19 |. Considering the long time behavior of the 
above expression, we have for x = 



2D \ 1/2 

m(0,t) 



1 + 



D 



One can also retrieve the spatial behavior with little effort. For small x we have 



(27) 



/ 2x 2 \ 1/2 

m(x,t)=m(0,t)+ — - +••• ,x<(Dt) 1/2 . (28) 



TlDt 



The large x behavior has two regimes: 



m(x,t) = l-(—z) exp -— +... , (Dt) 1/2 < z < At (29) 



, x (Xt\f2Dt\ 1 ' 2 f-x 2 \ x , , 

ro(M ) = l-^ H j(_) exp^j+... ,,»At. (30) 

It is interesting to note that apart from the natural diffusion scale (Dt) 1 ^ 2 , there is a larger 
'ballistic' scale At in the system, beyond which the disordering efficacy of the BA is much 
reduced, since it makes so few visits to these distant sites. There is no simple (i.e. single 
length) scaling form for m(x,t). 

Next we consider the continuum analogues of Eqs.(0)-([TI]). We define the average global 
magnetization (relative to its initial value) as 

oo 

M(t)= J dx [(0(x,O))-(0(x,t))] , (31) 

— oo 

which may be compared to the discrete version in Eq.(^). Integrating and averaging the 
continuum model (|I~5| ) yields 

^ = A(0(i?(t),t)), (32) 

which is to be compared with Eq.(||). This last equation indicates that we may explicitly 
find an expression for (<p(R(t) , £)) by calculating the time derivative of the spatial integral 
of m(x,t) = (<p(x,t)). This may be done at the level of the perturbation series (f2l|), from 
which one may show that the following relation holds exactly, for all times: 
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m(0,t) = {<j>(R{t),t)) , (33) 
which is the continuum analogue of Eq.(|7]). Finally, combining Eqs.(|3"2"D and (^) we have 

^ = An.(0,«). (34) 

Thus, the non-trivial relation between the rate of change of the global magnetization, and 
the mean of the magnetization density at the origin, is seen to be exact within the continuum 
model (which complements the exact relation fllCf ) found in the discrete framework). Similar 
results are easily obtained for all d G (0,2). In section V we shall derive a slightly more 
complicated form of (|3^) (involving the smearing function A/) which is appropriate for higher 
dimensional systems. Directly from Eqs. (|2"7|) and (^) we note that M{t) ~ \ft ; thus, the 
average number of corrupted bits in d = 1 increases as the square root of time. 

We now turn to spatial correlations in the system. These are most easily probed via the 
two-point correlation function 

C(x,t) = (<P(x,t)<p(0,t)) 

= (exp (-X J dt' [6{x - R{t')) + 6(R{1?))]j ^ , (35) 

where we have used the solution (|TjJ) in the second line. This average can be calculated 
using infinite order perturbation theory in A, just as was used to evaluate m(x, t). We write 

oo 

C(x,t) = jT i (-X) n c n (x,t) , (36) 

n=0 

with Co = 1. For n > 0, a given term c n can be explicitly evaluated by making integral 
representations of the n 5-functions, and performing the average over the paths R(t) (as 
described in Appendix B). Thus one has 

* n Tn ^ dki rdk 



c n (x, t)= I dn I dr 2 • • • / dr n / (1 + e~^ x ) ■ ■ ■ J (1 + e~^ x ) 





D " ' "' 

x exp ' 



z m=l \l=l 



2 



(37) 
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with the notation r n+ i = 0. This 2n-fold integral can be reduced using Laplace transform 
in time, such that the integrals over {/cj} may be performed, as described in Appendix C. 
The result is 



1 + exp 



(38) 



s{2sD) n / 2 

Summing over these functions with a weight of (— A) n and inverting the Laplace transform 
using a Bromwich integral 1 2011 , we have 

ds exp(st) 



C(x,t) 



(39) 



, 7 2tTI S 1 /2[l + S l/2 + eX p(- S l/2|x|)] ' 

where, as usual, the contour 7 is parallel to the imaginary axis, and to the right of any 
singularities. We have rescaled space and time as x = xX/D and t = t\ 2 /2D. 

This integral may be evaluated for large t in the following way. We re-express the integral 
as an expansion in powers of exp(— s 1 / 2 ^ |) (which is not the same as our original expansion 
in powers of A). So we have 



C(x,t) = ^(-l) n I n (x,i) , 



(40) 



n=0 



with 



ds 



3 st—ns 1 / 2 |a;| 



7 2vri sV2(l + s i/2) 



n+l 



,—ut 



f du ~~r77c 

27T J U 1 ' 2 




,—inv}-/ 2 \x\ 



1 /o 1 ~ 

„%nu ' \x 



+ 



(41) 



(1 + mV3)n (l_ m V2)r 

where the second line is the explicit form of the integral after integrating around the only 
singularity - a branch point located at s = 0. The integral over u may be simplified for 
t 1 to give 

1 x 



TTt 1 / 2 



{ixtyi' 



du 



It 



1/2 



COS 



W 2 



exp 



n 2 x 2 * 



(42) 
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In these rescaled units, we have from Eq.(p7|) m(0, t) ~ (nt) x l 2 (for large t). Thus, we may 
resum the functions {/„} to find 

C(x,t) ~ I n 2 i 2 \ 

scm * S (_1) exp ("in 

= ^ [l + ^(0, exp(-x 2 /2 J Dt))] , (43) 

where #4(5, z) is a Jacobi theta function (with norm q) Note, we have written the last 
line in unsealed variables, and we see that the ratio of the correlation function to m(0, t) 
does not depend on A for large times. The behavior of C (x, t) in the limits of large and small 
follows. For large x, the fields at the origin and at x will be uncorrelated, so that 
C(x, t) ~ m(0, t)m(x, t) ~ m(0,t), the latter result following since m(x,t) ~ 1 for x \/t. 
At the other extreme, as x — > 0, C(x,t) — > (0(0, t) 2 ). Referring to the exact solution of the 
continuum model, Eq. (|T6|) , one can see that the second moment of the magnetization density 
is actually given exactly by m(0, t;2A) (where the last (optional) argument indicates the 
parametric dependence on the spin-BA coupling). So for long times we take the expression 
for m(0,t;A) given in Eq.(PTp and replace A by 2A. Therefore (0(0, t) 2 ) ~ m(0, t;A)/2 for 
t 3> 1. Thus, the limits of the function C(x,t)/m(0,t) are 1/2 (for small x) and unity (for 
large x), which is naturally consistent with the analytic form given above in terms of the 
Jacobi theta function. In section VI we shall compare this expression with results from a 
numerical simulation of the discrete model described earlier. 

To complete our study of the properties of this system in one dimension, we shall consider 
the complete probability distribution P(0, x, t) of the magnetization density. We shall be 
able to calculate this exactly, since i) we can see from (|T6|) that the n th moment of the 
density is related to the mean density with a replacement A — >• n\; and ii) we have an 
exact expression for the mean density (albeit in the Laplace transform variable s). The first 
point is a fortuitous property of our continuum model which we should certainly exploit. 
The second property is less obvious. One might imagine that, given we know all density 
moments via the first, even the asymptotic form for the mean density would be sufficient 
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to calculate the probability distribution (for large times). This is not the case as we shall 
see - the complete analytic structure of m(x, t; A) is required in order to reconstruct the 
distribution V. 
We define V via 



r{<j>,x,t) = (s^-M^t))) 



(44) 



where <pR(x,t) is the stochastic field solution given in Eq. (|16[). We can re-express the 
^-function using a frequency integral, and expand in powers of the field as follows: 

dui 



— oo 

— oo n ~ u 

J 2tt ^ n\ 



(45) 



the last line following from the property i) alluded to above. 

So the Laplace transform (over time) of V is given in terms of the Laplace transform of 
m(x,t;nX). From Eq. (p5|) we have 



rh(x, s; nX) 



n\g(x, s) 
1 + nXg(o, s) 



<KM) -g(x,s) 



g{x,s J 



(46) 



sg(0,s) sg(0,s)(l + nXg(0,s)) ' 

where the second line follows from some algebraic manipulations. The first term is easily 
handled as it is independent of n. Thus the sum over n for this term (as is required in 
Eq.(f4~5"D) yields a factor of exp(iu;) which finally yields a factor of 5(1 — 0) when integrated 
over io. The second term is more interesting. Details of how to perform the sum over n and 
the frequency integral may be found in Appendix D. The final result for V(4>,x, s) reads 



v 9(0, s) - g{x, s) g(x,s) 1 

V{(p,X,S) = — — r 6(1 — (f>) + TTTT-x?— TT ex P 



sg(0, s) 



^(O, S ) 2 sA0 



Xg(0,s) 



lo S It 



(47) 
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This form for the Laplace transform of the probability distribution may be easily generalized 
for any d G (0, 2). Finally we must invert the Laplace transform. To this end we note that 

1 



g(x,s) 



cxp 



d) |x| 



Inserting ((48]) into ( 
V(<t>,x,t) = 6(l- 



(2Ds) 1 /2 

and inverting the transform we have our final result 



(48) 



erf 



\x\ 



1 



1 



exp 



\x\ log(0) 

(2Dty/ 2 ~ 2AtV2 



(49) 



where erf (z) is the error function jBJ, and we have defined A = A/ (2D) 1 1 2 . This is illustrated 
in Fig. 2 for x — 1 (and D = 1, A = \/2), and three different times corresponding to 
x 2 > 2Dt, x 2 ~ 2Dt and x 2 < 2£>t. 

In particular, the probability distribution for the magnetization density at the origin 
takes the form 

(log(0)) 2 " 



exp 



4A 2 t 



(50) 



(7rt)V2 A 

which is a pure log-normal distribution. This result is very revealing, as it shows that the 
fluctuations in this system are extreme. For instance, we have already seen that the mean 
value of the density (at the origin) decays as m(0, t) ~ However, if one asks how the 

typical (or most likely) value decays, one can see from (|50| ) that (0(O,t)) typ ~ exp(— X 2 t/2). 
Thus, as time proceeds, the typical value of decays to zero exponentially fast, whilst the 
mean decays slowly as 1/ \ft. This is possible because the log-normal distribution has a long 
tail, extending out to the extreme value of <fi = 1. In fact the end point of the distribution 
(i.e. V(l, 0, t)) also decays as 1/y/t which is consistent with the known persistence properties 
of a random walker in d = 1 (namely, the probability of a walker never having returned to the 
origin after time t decays as 1/y/i). In Fig. 3 we illustrate V((fi, 0, t) for three different times. 
As a final remark, we note that if we erroneously use the asymptotic form ( |27|) for m(0, t; nX) 
to build the distribution function, we find that d$P(<\), 0, t) is equal to d^,8(4>) — 5(1 — 
thus emphasizing the fact that we need the entire analytic form of m(0, t; nX) to successfully 
construct the distribution V . 
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This ends a rather long section on the analytic properties of the continuum model in 
d — 1. In the next section we shall briefly study the case of higher dimensions, and then in 
section VI we shall compare our results with some numerical simulations performed on the 
original discrete model. 



V. RESULTS IN HIGHER DIMENSIONS 

As we mentioned several times in the previous two sections, the continuum theory re- 
quires some regularization for d > 2. This can be most easily (and physically) accomplished 
by smearing the ^-function interaction between the BA and the spins. Thus, our continuum 
model takes the form 

«9 t 0(r, t) = -A0(r, *)A,(r - R(t)) , (51) 

where Aj(r) is a normalized function which is sharply peaked over a region of linear dimension 
I around the point r. [A good choice would be A; ~ l~ d exp(— r 2 // 2 ).] 

In this section we shall analyze some basic properties of Eq.(|5TD for general Aj. Then 
we will use a somewhat more crude approach to estimate the mean magnetization density 
(at the origin) function of time in d = 2. 

Although we have generalized our continuum theory somewhat, we can still make sub- 
stantial headway by first integrating the equation of motion, and then using infinite order 
perturbation theory. The first step yields 



(r, t) = exp 



-A J dt' A,(r-R(0) 



o 



(52) 



whilst the second consists of expanding this equation in powers of A and averaging term by 
term: 



m 



r,t) = (0(r,t)) = ]T(-ArxnM) , (53) 



n=0 



where Xo( r , t) = 1, and for n > 
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Xn(r,t) 



nl 



J dr A,(r-R(r)) 



(54) 



By making a Fourier representation of the interaction function A;, performing the average 
over paths (see Appendix B), and finally Laplace transforming in time, we arrive at 



Xn{r,s) - < ll l , o _ w j- f exp[ZK n rj 



^t X J (27T) 



(55) 



with the convention k = 0. 

We shall be concerned with two quantities. The first is a smeared mean magnetization 
density near the origin, and the second is the global magnetization. These are given by 



m 8 (0,t) = J d d r A z (r)m(r,t) 



(56) 



and 



M{t)= d d r [<0(r,O))-<0(r,t)>] 



(57) 



respectively. We shall prove that 



dM{t) 
dt 



Xm s (0,t) , 



(58) 



which is the smeared analogue of the global versus (strictly) local relation ([34) we proved 
for d < 2. Comparing Eqs. (p3|) and (^) we see that 



n=0 



m s (0,t) = £(-A)X(0 
where in Laplace transform space 



(59) 



Us) = Jd d r M^Ms) = - ( ft / ' 

lm=l 



<i /c m A;(k m k m _i^ 

(27p (^ + fcS 



A ; (k, 



(60) 



Similarly we have 



dM(t) 
dt 



E("A) n Cn(t) 



(61) 



n=0 



with 
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Cn(t) = J t J d d r X n(r,t) . (62) 
For n > 0, Xn( r , 0) = 0, so that 

Us) = s fd d rxn(r,s) . (63) 



Using Eq.(^) we may evaluate the above integral to give 

d k m A;(k m k m _ 



us) H n / ^ "T + ~*sr nM • (64) 



vm=l 

Performing the integral over k n immediately yields 



Cn(«) = ^n-i(s) • (65) 



Thus comparing Eqs. (|59|) , (|6l|) and (|65|) we see the validity of relation Eq.([58|). 

As a corollary, by integrating the averaged equation of motion floTf ) over space, we have 

= A / d d r (Mr - R(t))0(r, t)) , (66) 

When compared with Eq.(j58|), the above relation gives us 



m s {t) = J d d r A,(r)(0(r,t)> = f d d r (A,(r - R(t))0(r, t)) , (67) 

which is the smeared version of the local relation ( p3|) proven in section III. We note that, 
although we have been concerned with a sharply peaked interaction function, the relations 
(EB|) and (K7p hold for an?/ function A(r). 

This ends the more rigorous part of the present section. In the remainder we shall just 
mention some explicit results for the mean local magnetization density (at the origin), which 
are obtained with a cruder regularization. 

The difficulty with making headway using the smoothing function, is that the n-fold 
integrals over the A/'s are intractable (unless one can find a particularly 'friendly' form 
for Aj.) As an alternative approach, we return to the sharp Dirac 5-function as used in 
section IV. We remarked that the n-fold convolution integrals were divergent due to the 
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non-integrability of g(0,t) for d > 2. To evade this difficulty we can simply impose a cut- 
off into the integration limits. This is closely connected to introducing a microscopic time 
scale into the temporal correlations of the BA. Such an regularization procedure was used in 
Ref . |J , and the results so obtained were shown to be equivalent to previously known exact 
results 0. So we shall use the same procedure here, but with due caution. 

First, we consider d = 2. In a precisely analogous way to the calculation in section IV, 
we expand the field solution fll6|) in powers of A and average term by term. Thus, we have 
(cf. Eqs. (0) - (H) 

<x> 

m(0, t) = (0(0, *)> = £ (-X) n Xn(0, t) , (68) 

n=0 

where Xo(0, t) = 1, and for n > 

t— to Tl—to T"n—1— *0 

Xn (0,t) = J dn J dr 2 --- J dr n g(Q, t x - r 2 ) • • • g(0, r n _i - r n )g(0, r n ) , (69) 

to to to 

In this case, the probability distribution of the BA at the origin has the form g(0,t) = 
(2'irDt)^ 1 . Note we have inserted the microscopic time regulator to in the limits of the time 
integrals. Our strategy is to evaluate the time integrals one by one, keeping only the most 
singular term at each step. We shall use the general result (for t ^ t ) 



7° dr [log(T-Ao)] m „ (m + 2) [\og(t/t )r^ 

J (t-r)r ~ (m+ 1) t ' { ' 



(t-r)r (m+1) 

to 

Therefore we have the dominant contribution 



^(°'*)~ (2^)" ' ^ 



Inserting this result into Eq. (R3) and summing over n we have the asymptotic form 



m(0,t) 



A / 1 

1+ M log U 



(72) 



Thus the magnetization at the origin does decay to zero for large times, but logarithmically 
slowly. From the relation fl58|) we see that the mean global magnetization increases as 
M(t) ~ t/\og(t). 
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The same kind of analysis can be repeated for d > 2, and one finds that m(0, t) saturates 
to a constant for large times, which implies that M(t) ~ t for large times. These results are 
easily understood from the recurrent properties of the BA (i.e. a random walker returns 
to its starting point with probability one, only for d < 2). It would be more interesting to 
derive the distribution of the magnetization density for d > 2, but this requires the more 
careful regularization method involving the smoothing function Aj(r) and thus lies beyond 
the scope of the present work. 

[In d = 2 one can capitalize on the slightly crude result ( fT2|) obtained using the cut- 
off to, combined with the exact property (0(0 , t; A) n ) = m(0, t;n\) to derive a form for 
the distribution function V((j),0,t). Such an approach yields V((f),Q,t) = /3(t)0^ -1 , with 
(3 = 2nD/\ \og(t/t ). However, this result is not to be taken seriously, since we need the 
whole analytic structure of m(0,t;n\) in order to derive P, as evidenced in the previous 
section.] 

VI. NUMERICAL SIMULATION 

We have performed extensive numerical simulations of the discrete model, as defined in 
section II, in order to test the results obtained in the last two sections from the continuum 
theory. In all of the simulations for which we present results, we have set the hopping rate 
p of the BA, along with the flipping probability q, to unity. We have experimented with 
decreasing the flipping probability, and have found that its only effect is to renormalize the 
effective spin-BA coupling A, such that A oc q, as expected. 

Most of our results are obtained from a one-dimensional chain of sites. The chain length 
is unimportant, so long as one ensures that the BA has never touched the edges in any of 
its realizations up to the latest time at which data is extracted. Generally we average over 
between 10 6 and 10 8 realizations (or runs) depending on the desired quality of the data. 
Such simulations required a few days on a DEC Alpha 233 MHz workstation. In a given 
run, at each time step the BA is moved left or right with equal probability and the spin 
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it leaves behind is nipped. Each run starts with the same initial configuration; namely all 
spins up, except the spin at the origin (which is the starting site of the BA) which is pointed 
down. (This means that all spins are up after one time step, since the BA has moved away 
from the origin and flipped the down spin.) 

In Fig. 4 we show the measured mean magnetization density at the origin, along with the 
quantity J2r ©(A t, R) (cf. Eq.(@)). They are seen to be identical thus confirming relations 
Eq.(0) and its continuum counterpart fl33[ ). The solid line is the asymptotic prediction fl2~7[) 
from the continuum theory. It is seen to be in good agreement with the data, as the line 
has a slope of (-1/2). From the fit of this log-log plot we can read off the effective value of 
A, since from fl2~7|) the prefactor of l/y/i is given by (2/7rA 2 ) 1//2 . (The diffusion constant for 
the lattice random walk is unity.) We have fitted the data to c/y/i with c = 0.40(1), which 
yields A = 1.99(1). In Fig.5 we plot the small xj \fi dependence of m(x, t) on a log-log scale. 
The data is well fitted by the prediction given in Eq. (p8|) . In Fig. 6 we plot 1 — m(x, t) versus 
x/yt. Note that good data collapse is found for intermediate values of xj \ft. We have been 
unable to numerically probe the ballistic scale Xt. (Note, the theoretical curves shown in 
the last two figures are plotted with no free parameters.) 

In Fig. 7 we plot the discrete time derivative of the total number of down spins N(t) 
(which is M(t)/2), along with m(0,t). The two curves are indistinguishable within the 
numerical noise, thus confirming the global/local relation fllPf ). This also provides secondary 
confirmation of the continuum form of this relation Eq.fl34|) with A ~ 2. 

In Fig. 8 we plot the ratio of the measured two point correlation function (divided by 
m(0,t)) versus x/yt. Note that it varies from 1/2 (at small x) to unity (at large x) as 
expected. The data from two different times is shown, and one sees excellent agreement with 
the theoretical prediction (f43|), which is plotted with no free parameters. This agreement 
provides very strong evidence for the validity of our whole continuum approach. 

Briefly, we mention simulations in d — 2. Higher dimensional simulations are not too 
difficult as one is only ever moving the single BA at each time step. In Fig. 9, we show 
mo(t) and J2r ©(R-> ^ R-)- The data for the two functions are identical verifying the discrete 
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relation (0) in two dimensions, as well as confirming the continuum result (|67|). We have 
plotted the inverse of these functions against log(t) in order to compare with the theoretical 
prediction ([T2|). Again, good agreement is found, thereby confirming the less rigorous method 
by which the two dimensional result was obtained. 

Finally we mention our attempt to measure the probability distribution of the coarse 
grained magnetization density at the origin (in d = 1), which was found from the continuum 
theory to be a log-normal distribution. Clearly it does not make sense to measure moments 
of the spin at the origin, since the odd (even) moments are equal to mo(t) (unity). Therefore, 
we define a coarse grained magnetization over a patch of spins. If the patch is taken too 
small, the coarse-graining will be ineffective, whilst if the patch is taken too large, the BA 
will take a long time to leave the patch, and the asymptotic behavior will be numerically 
inaccessible. So we have compromised and have used a patch containing 21 spins. We have 
binned the patch magnetization from 10 8 independent runs and generated the histograms 
shown in Fig. 10. Note, that because the patch size is modest, the histograms have non-zero 
weight in the negative region, in contrast to the strict continuum limit. However, we do 
see that for near unity, the histograms have a robust tail, which is the signature that 
extreme fluctuations are important. 

VII. CONCLUSIONS 

In this paper we have introduced and analyzed a simple model of data corruption due 
to a Brownian agent (BA). In section II we introduced a discrete version of the model, 
which consists of a BA flipping bits (or spins) on a lattice. The model is non-trivial since 
the value of a given spin depends very sensitively on the path of the BA (i.e. whether the 
spin has been visited an odd or even number of times). We presented a master equation 
formulation of the model and derived an equation of motion for the conditional average 
of the magnetization density. In the continuum limit, this quantity was seen to satisfy an 
imaginary-time Schrodinger equation (ITSE) for a particle in a repulsive 5-function potential. 
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Higher-order conditional averages also satisfy ITSE's with an additional repulsive ^-function 
potential for each spin being averaged. We also proved that there is an exact proportionality 
between the rate of change of the mean global magnetization M(t) (as defined in Eq.(|8|)) 
and the mean magnetization density at the origin, thus revealing a non-trivial statistical 
relation between a global and a local quantity. 

In section III we recast the discrete model as a stochastic cellular automaton (SCA). 
From the local rules of the SCA, we postulated a particularly simple continuum theory 
written in terms of a coarse-grained magnetization density 0(r, t) (which is a functional of 
the walk R(t) of the BA). We showed that averages of this stochastic field may be recast 
as integrals over the final BA position, the integrands of which are imaginary-time path 
integrals (ITPI's). These ITPI's were shown to describe quantum systems of repulsive 5- 
function potentials identical to those revealed through the ITSE formulation of the master 
equation. In this way we have confirmed that the continuum theory is a good representation 
of the underlying discrete model. [It is worth mentioning that averages of 1/0 correspond 
to the quantum mechanics of an attractive 5-function potential. In this case there is the 
possibility of a sharp transition in behavior for d > 2 as one varies A, since bound states 
only exist for a sufficiently attractive well.] 

In section IV we thoroughly examined the properties of the continuum theory for d = 1. 
First, we derived an exact expression for the evolution of the magnetization density m(x, t), 
and found its spatial variation for small and large x. In particular we found m(0, t) ~ 1/ X\/i 
for large times; and also the existence of both a diffusive and a ballistic scale in m(x,t). 
Using the exact result for m(x, t) enabled us to prove that the rate of change of M(t) is 
proportional to m(0, t) for all times (with proportionality constant A) which is the continuum 
analog of the global/local relation that was proved in section II. We then studied the two- 
point correlation function C(x, t). An exact expression was found for the Laplace transform 
of this quantity, from which we were able to extract its long-time behavior, which is expressed 
in terms of a Jacobi theta function as shown in Eq.([J3]). The correlation function has an 
asymptotic scaling form; namely the ratio C(x, t)/m(0, t) depends only on x/\/t. Finally in 
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section IV, we examined the probability density V of the magnetization density. Using a 
fortuitous property of the original continuum theory which enables us to represent the n th 
moment of the magnetization density in terms of m(x,t), along with our exact expression 
for this latter quantity, we were able to derive an exact expression for V(<f), x, t) as shown in 
Eq.(^Uj). In particular, setting x = reveals that V(<f), 0, t) is a pure log-normal distribution. 
This last result emphasizes the importance of extreme fluctuations in this system. For 
instance, m(0,t) decays as l/Vt, whilst the typical value of the magnetization density at 
the origin (i.e. the mode of V) decays as exp(— X 2 t/AD). 

In section V we briefly examined higher dimensions. For d > 2 it is necessary to regularize 
the continuum theory, most appropriately with a smeared interaction between the BA and 
the spins. Using an arbitrary interaction function A;(r) (which has a linear scale /) we were 
able to prove a more general form of the global/local relation, as shown in Eq.fl58|). We then 
concentrated on d = 2, and using a cruder regularization (namely, introducing a microscopic 
correlation time to) we were able to derive an expression for the asymptotic decay of m(0, t). 
This decay (cf. Eq. ([72]) ) is seen to be logarithmically slow. 

In section VI we presented our results from extensive numerical simulations of the origi- 
nal lattice model. Most of our numerical work examines the case of d — 1. We measured the 
spatial and temporal variation of the mean magnetization density, the mean global magne- 
tization, and the two-point correlation function C x (t). In all cases we found excellent agree- 
ment between our data and the theoretical predictions arising from the continuum model. 
In particular, the agreement between the theoretical form for C(x, t) and the numerical data 
is very satisfying, as there are no free parameters to adjust. We also measured m Q (t) in two 
dimensional simulations, and found good agreement with the predicted logarithmic decay. 
Finally we attempted to measure the probability distribution of the magnetization density, 
by binning the magnetization of a spin patch (containing 21 spins) from 10 8 realizations. 
The results are not of high enough quality to directly compare with the derived log-normal 
distribution; however, we do see clear evidence of a long robust tail of the histogram for 
values of <fi near unity, which is a clear signature of the importance of extreme fluctuations. 
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In conclusion, we have introduced and solved a model in which a BA interacts with a 
bimodal environment (i.e. a medium containing two types of particles, spins, bits, etc.). 
Our primary application has been an environment composed of bits of data, which the BA 
steadily corrupts. We have been interested in the statistical correlations and fluctuations 
of the disordered medium, and our exact results (arising from a simple continuum model) 
reveal the correlations to be non-trivial and the fluctuations to be extreme in nature. These 
statements are made quantitative by the form of the two-point correlation function, and 
the probability distribution of the density of corrupted bits. There are many directions 
for future work, foremost amongst which are: i) calculating the distribution V in d = 2 
and determining its sensitivity to the smearing function A;, ii) investigating autocorrelation 
effects in one dimension, iii) refining the model to take into account relevant factors (like 
a back-coupling between the environment and the BA, or an asymmetry in the flipping 
probability) which will arise if one tries to make a stronger connection to real processes. 

The model may also be seen to be a very simplified version of other systems. For instance 
the bistable medium can be taken to be composed of two chemical species A and B (with 
vanishingly low mobility) and the BA to be a high mobility catalyst, inducing a reversible 
reaction between A and B (and vice versa). Alternatively we can think of the BA as a 
wandering impurity in an ionic crystal (such as an anion or cation vacancy in NaCl) or a 
semiconductor compound (such as Zn in GaAs), which has a small probability of reordering 
the local bi-atomic structure as it passes through a given unit cell [^1] . 



The data corruption process appears to us the most interesting application, as well 
has being the most potentially relevant. This is especially true given the enormous efforts 
dedicated to creating memory storage devices of ever-decreasing size. Such miniaturization 
will lead to new causes of soft error production [|12|], amongst which will inevitably be found 
the Brownian agent. 
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APPENDIX A: 



In this appendix we outline the solution to the discrete equation (0) for the conditional 
average 0(r, t | R). This is achieved by use of discrete Fourier and Laplace transforms 
defined via 

oo 

F k \ R {£ z \Mv,t | R)]] = £0(r,n<ft | R)exp(zk-R) , (Al) 

n=0 R 

(with t = n5t). By self-consistently determining the function 0(r,t | r), one can explicitly 
solve for the double transform of 0. Inverting the discrete Fourier transform and summing 
over R one finds (with the particular choice of initial condition (^|)) 

r 1 - 2(1 - z) J d d k E(k, z) 1 
£ £ "' [e(R ' ' 1 R)1 = (1^) { 1 + f J #k J(k, ,) Ei e--- ) ■ (A2) 

where J BZ indicates integration over the first Briulloin zone, and 

E(K z) = 1 1 - z + | B 1 - exp(-zk ■ 1)] | . (A3) 
Similarly one can find an explicit expression for Z)£ z it[0(O, t | R)]. The expressions can be 

R 

shown to be equal thus proving the relation (0) as given in the main text. 



APPENDIX B: 



In this appendix we outine the procedure for averaging the expression in Eq.(^). First, 
we time order the integrals, thereby absorbing the factor of l/n\. Then, making an integral 
representation of each Dirac 5-function, we have 



Xn(x,t) = J dr x j dr 2 ■ ■ 



x (exp [— ikiR(ri) 



dr n I dhe iklX 
■ ■ ■ — ik n R(r n 

Using the solution of Eq. flT3|) , we may rewrite the average in ( |B1| ) as 

i(h + --- + k n ) J dr£( 
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The Gaussian average over the noise may now be performed, and the above expression 
reduces to 

exp j-^ [kUn - r 2 ) + (h + k 2 ) 2 (r 2 - r 3 ) + • • • + (h + ■ ■ ■ + k n ) 2 r n ]\ • (B3) 

We now insert this expression back into Eq . flBl] ) , and make the change of variables k[ = 
k\, k' 2 — k\ + k 2 , ■ • • i k' n — ki + ■ ■ ■ + k n . The integrals over {k'j} are easily performed and 
we arrive at Eq.(^) in the main text. 



APPENDIX C: 



In this appendix we outline the evaluation of the 2n-fold integral for c n (x,t) given in 
Eq. (p7|) . We note first that the time integrals have the form of a ra-fold convolution. Thus, 
we can Laplace transform the expression to find 



dki 
27 



1 + e- lklX ) 



x n i[s+(jD/2)(fci+ ... +U2] 



(Cl) 



We make the change of variable k[ = ki, k' 2 = k\ + k 2 , ■ ■ ■ , k' n = ki + • • • + k n , and rearrange 
the integrals to give 



Crj I X . S J 



i/f [s+ p /2 )^r'-/f 1.HD/2KI- 1 

n 

x [] [l + expH(C - Ki-i)z]] , 



m=l 



(C2) 



with the understanding that k' = 0. We now multiply out the product which gives us 2" 
terms. These are grouped into n + 1 sets, the m th set containing terms which are equal 
after integration. Thus each term in the m th set (within the integrals) may be taken to 
contain m factors of e tkjX (where j = 1, • • • ,m) and (n — m) factors of unity. Using the 
integral 



dk 



Akx 
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we can reduce Eg. (P2|) to the form 

1 n 

b m=0 



Performing the binomial sum gives Eq.(|38|) in the main text. 
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(C4) 



APPENDIX D: 

In this appendix we outline the derivation of Eq . (f47|) from Eqs.([15|) and fl46|) . The only 
non-trivial aspect of the derivation is the frequency integral and sum (over n) of the second 
term in Eq . fl4"6|) . Ignoring the prefactor of that term (namely g(x, s)/ sg(0, s)) we must 
evaluate a quantity Q(s, <p), which has the form 

" du ™ (iu) n 1 



QM) ^ / E 



- n! (1 + nA£(0, s)) ' 
In order to perform the sum, we introduce the integral representation 

no 

1 



{l + n\g(0,s)) 



du exp — u(l + n\g(0, s)) 



A 

due' u ( exp (iue~ uX§ 

J 2tc v 

-oo 



The sum over n now reconstitutes an exponential function, and we have 
The integral over u is easily done to give 

oo 

Q(s,<j>) I da < " 8 (o -"^ 



Finally, changing variables to v = e uX ^ we have 



Q(s, 



\g(O,s)0 



exp 



\g(0,s] 



log 



(Dl) 



(D2) 



(D3) 



(D4) 



(D5) 



Using this result for Q in conjunction with Eqs. (|45|) and (|46D we have ( f47[) in the main text 
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FIGURES 
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FIG. 1. Illustration of the data corruption process for d = 2. The initial uncorrupted state is 
shown on the left, with the BA represented by the filled circle. On the right we show a typical 
walk of ~ 20 steps. The BA flips a spin with each visit, so those spins visited an even number of 
times are restored to their original value. 
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FIG. 2. V((j),x,t) versus <fi, as given in Eq.(|49|), with x = 1,A = \/2, for three different times, 
i = 0.5, 1.0, 1.5. The thick vertical line represents the 5- function at <f> = 1. 
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FIG. 4. Log-log plot of mo(t) (diamonds) and ^2 R @(R,t, R) (plusses) versus time from d = 1 
numerical simulation. The solid line is the asymptotic theoretical prediction (pTj). 
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FIG. 7. Log-log plot of the time derivative of the total number of minus spins (diamonds) and 
of TOo(t) (plusses) versus time in d = 1. The solid line is the theoretical prediction (p7[). 
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